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An XMM-Newton observation of the Narrow Line Seyfert 1 
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ABSTRACT 

XMM-Newton observations of the NLS1 Markarian 896 are presented. Over the 2- 
10 keV band, an iron emission line, close to 6.4 keV, is seen. The line is just resolved 
and has an equivalent width of ~ 170 eV. The broad-band spectrum is well modelled 
by a power law slope of T ~ 2.03, together with two blackbody components to fit the 
soft X-ray excess. Using a more physical two-temperature Comptonisation model, a 
good fit is obtained for an input photon distribution of kT ~ 60 eV and Comptonising 
electron temperatures of ~ 0.3 and 200 keV. The soft excess cannot be explained purely 
through the reprocessing of a hard X-ray continuum by an ionised disc reflector. 
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1 INTRODUCTION 

Narrow Line Seyfert 1 galaxies (NLSls) are a subset of 
AGN with particularly narrow Balmer lines, these being only 
slightly broader (H/3 < 2000 km s _1 ; Osterbrock & Pogge 
1985) than the forbidden lines; in this respect they are sim- 
ilar to Seyfert 2s. However, NLSls also show an [Oiii] to 
H/3 ratio of >3, which is one of the distinguishing features 
of Seyfert 1, rather than Seyfert 2, type galaxies (Shuder & 
Osterbrock 1981). The third characteristic is the frequent 
presence of strong Fe II (or higher ionisation iron) emission 
lines; these, again, are more usually seen in Seyfert Is than 
Seyfert 2s. 

Observations with ROSAT found that NLSls tend to 
have steep 0.1-2.4 keV X-ray slopes, leading to soft luminosi- 
ties which are higher than can be explained by the reprocess- 
ing of the hard X-ray continuum. They also have increased 
variability in this region, with doubling timescales as short 
as a few minutes (e.g., Boiler, Brandt & Fink 1996). A pos- 
sible explanation for the extreme properties of these AGN is 
that they contain relatively low mass black holes, which are 
accreting at a high (near-Eddington) rate (Pounds, Done & 
Osborne 1995). 

In this paper, we present XMM-Newton observa- 
tions of Mrk 896 (also known as MCG-01-53-008 and 
IRAS F20437-0259). With an H/3 FWHM of 1135 km s" 1 
(Veron-Cetty, Veron & Gongalves 2001), Mrk 896 lies to- 
wards the centre of the range of NLSls line-widths (~ 500- 
2000 km s" 1 ) and is a low redshift (z = 0.0264) object, with 
a 5GHz flux density of 38 mjy (Bicay et aJ. 1995). Boiler et 
al. (1996) measure T to be 2.6 ± 0.1 over the ROSAT band 



of 0.1-2.4 keV, while it was found to be 2.82^°;" during an 
earlier observation (Boiler et al. 1992). 



2 XMM-NEWTON OBSERVATIONS 

Mrk 896 was observed twice in revolution 355 (2001 Novem- 
ber 15), as part of a sample of NLSls. All instruments (EPIC 

- European Photon Imaging Camera - Striider et al. 2001, 
Turner et al. 2001; RGS - Reflection Grating Spectrometer 

- den Herder et al. 2001; OM - Optical Monitor - Mason et 
al. 2001) were used over the second observation. During the 
first observation, the EPIC instruments were closed; data 
were obtained for the RGS and OM, however. The total 
exposure times for the various instruments were as follows: 
MOS - ~ 10 ks; PN - ~ 8 ks; RGS - ~ 20 ks; OM - ~ 16 kS 

The pipeline-produced EPIC event-lists were filtered us- 
ing xmmselect within the SAS (Science Analysis Software); 
single- and double-pixel events (patterns 0-4) were chosen 
for the PN data, while the range 0-12 was used for the MOS 
cameras. Source spectra were extracted from the images us- 
ing a circular region of 45 arcsec radius. This same region 
was then moved to an adjacent area of 'blank sky', to obtain 
a background spectrum. 

The RGS 1 spectra for the two observations were co- 
added, as were the RGS 2 datasets. The spectra were 
binned, using the FTOOL command GRPPHA, to provide a 
minimum of 20 counts per bin for the EPIC spectra. The 
Xspec vll.0.1 software package was then used to analyse 
the background-subtracted EPIC and RGS spectra, using 
the most recent response matrices, (rgsrmfgen was run 
within the SAS to obtain the relevant RGS responses.) 
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Throughout this letter, Ho = 50 km s Mpc and 
qo = 0.5 are assumed; unless stated otherwise, errors are 
given at the 1 a level. 



3 SPECTRAL ANALYSIS 
3.1 EPIC data 

As is conventional, the first model fitted to the joint 
MOS+PN data was a single absorbed power law (Nh fixed 
to the Galactic value of 3.89 x 10 20 cm" 2 ; Murphy et al. 
1996), over the entire 0.2-10 keV bandpass. This provides a 
poor fit, with xt = 1168/720, mainly due to a strong upward 
curvature below ~ 2 keV. Constraining the model solely to 
the 2-10 keV band, however, produces an excellent fit for a 
photon index T ~ 1.96 (\l = 261/266; fit 1 in Table[TJ. 

There were seen to be positive residuals to the fit around 
6.4 keV, in both the MOS and PN data, so a Gaussian emis- 
sion component was added. Allowing both the energy and 
width of the line to go free, the best fit was obtained for a rest 
energy of (6.36 ± 0.06) keV, with a = (0.124 ± 0.075) keV. 
The equivalent width was measured to be 182 eV (fit 2 in 
Table 0. This improved the fit by A\ 2 of 13 for 3 degrees 
of freedom, corresponding to a probability of >99 per cent. 
A second Gaussian did not improve the fit significantly. It 
should be noted that freezing the energy to 6.4 keV, for neu- 
tral iron emission, gave an equally acceptable fit. Also, if the 
line is fixed to be unresolved (a — 0.01 keV), the fit is only 
worse by A\ 2 of one, for one degree of freedom; this line 
has an equivalent width of ~ 140 eV. This is consistent with 
lines seen in other objects (O'Brien et al. 2001; Reeves et 
al. 2001; Pounds et al. 2001; Kaspi et al. 2001; Lubinski & 
Zdziarski 2001). However, if we assume there is no intrinsic 
narrow line, a diskline model (Fabian et al. 1989) can be 
tried (fit 3); this was found to be a good fit, with an incli- 
nation angle of 24 ± 15° for the disc (inner radius fixed at 
a typical value of three Schwarzschild radii; emissivity in- 
dex = —2). In this case, the narrow component corresponds 
to the blue 'horn' of the profile. It is, unfortunately, not fea- 
sible to decide with statistical significance between between 
the Gaussian and diskline models. 

Extrapolating the 2-10 keV power law down to 0.2 keV 
revealed an obvious broad soft excess (Fig. 0. To model 
the broad-band spectrum, the iron line width and energy 
were fixed to the 2-10 keV values and blackbody (BB) com- 
ponents used to parametrize the soft excess. It was found 
that the breadth of the excess required two BBs to fit the 
observed spectrum (fit 5 in Table 0. 

Since the previous ROSATpwpers had modelled the soft 
excess region using a power law, it was decided to try fitting 
the XMM soft X-ray spectrum with a simple power law. Over 
the 0.1-2.4 keV ROS AT h&nA-pass we get V = 2.55 ± 0.01, in 
close agreement with the value given by Boiler et al. (1996). 
Using either two separate power laws or a broken power 
law model over the 0.2-10 keV bandpass (fits 6 and 7 in 
Table0 provides statistically acceptable fits, although they 
are much worse than for the blackbody parametrisation (or 
thCompfe; see below). Thus the soft X-ray spectrum shows 
intrinsic curvature. 

The broad-band best-fit model (fit 4) gives a 0.2-10 keV 
flux of 7.85 xlO -12 erg cm -2 s _1 , which corresponds to a 
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Figure 1. The 0.2-10 keV MOS (black) and PN (grey) spectra 
of Mrk 896. The 2-10 keV power law has been extrapolated down 
to lower energies, showing the presence of a strong soft excess. 
Positive residuals can be seen just above 6 keV, indicating an 
emission line. 



luminosity of 3.35 xlO 43 erg s _1 ; ~ 30 per cent of this lu- 
minosity is due to the multiple blackbody components used 
to model the soft excess. The 2-10 keV luminosity was mea- 
sured to be 9.67 xlO 42 erg s" 1 . Over the ^OSAT band, the 
power law model gave a luminosity very similar to the value 
given in Boiler et al. (1996) but about half that seen in an 
earlier R OS AT observation (Boiler et al. 1992). 

3.1.1 Variability Analysis 

Figure |21 shows the MOS and PN 0.2-10 keV lightcurves for 
Mrk 896 during the XMM-Newton obervation. There is a 
±20-25 per cent variation in count rate over a few thousand 
seconds, which is consistent with other NLSls (e.g., Boiler, 
Brandt & Fink 1996). Boiler et al. (1996) saw a variation 
by a factor of 1.9 in just under 300 ks; here, we find a much 
faster variation, of ~ 1.5 in 3 ks, although even this more 
rapid change only corresponds to a radiative efficiency value 
of ~ 0.3 per cent (Fabian 1979). Over this short observation, 
we find no evidence for a difference in the variability ampli- 
tude when comparing the 0.2-1 keV and 1-10 keV bands; the 
calculated hardness ratio is fully consistent with being con- 
stant throughout the observation (lower panel of Figure 
Using a cross-correlation function, no significant time delay 
between the hard and soft bands is found, with an upper 
limit of < 270 seconds. Assuming a black hole mass of 10 6 
M© (see section l3.1.2H . this sets a limit of < 27 Schwarzschild 
radii for the size of the emitting region. 

3.1.2 Comptonisation fits 

Although modelling the spectrum with a power law and BB 
components produces a good fit, it has no physical basis. 
Using the relationship between the temperature of the ac- 
cretion disc and the black hole mass (e.g., Peterson 1997) 
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Table 1. 

Fits to the XMM-Newton data. 

a Rest-frame energy of the emission line. 6 Intrinsic line-width. c Blackbody temperature. d Break Energy of 1.08 ± 0.04 keV. * frozen. 



Fit 


Range 
(keV) 


Model 


r 




(keV) 


(keV) 


EW 

(eV) 


kT c 
(keV) 


X 2 /dof 


1 


2-10 


PL 


1.96 


± 0.06 










261/266 


2 


2-10 


PL+GA 


2.04 


± 0.04 


6.36 ± 0.06 


0.124 ± 0.075 


182 ± 65 




246/263 


3 


2-10 


PL+DISKLINE 


2.06 


± 0.04 


6.42 ± 0.12 




244 + 40 
Z44 -115 




245/263 


4 


0.2-10 


PL+GA+BB 


2.16 


± 0.02 


6.36-f 


0.124^ 


218 ± 50 


0.088 ± 0.002 


713/717 


5 


0.2-10 


PL+GA+2BB 


2.03 


± 0.04 


6.36^ 


0.124^ 


165 ± 48 


0.084 ± 0.003 
0.224 ± 0.022 


698/715 


6 


0.2-10 


BKNPL+GA 


2.65 
2.15 


± 0.02 d 
± 0.02 


6.36-f 


0.124/ 


222 ± 54 




740/717 


7 


0.2-10 


2PL+GA 


2.94 
1.71 


± 0.09 
± 0.10 


6.36-f 


0.124/ 


152 ± 48 




741/717 
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Figure 2. Lightcurves showing the 0.2—10 keV continuum vari- 
ability of Mrk 896 over the 10 ks observation, using 125 second 
bins. The PN camera was switched on approximately 2400 sec- 
onds after the MOS instruments. The lower panel shows the hard- 
ness ratio between the 0.2—1 and 1-10 keV bands. 



(where M is the mass accretion rate, MEdd is the Ed- 
dington accretion rate, Ms signifies the mass of the central 
black hole in units of 10 s Mq and R sc h is the Schwarzschild 
radius), the temperature expected for the inner-most part 
of the accretion disc (r = 3 K sc h), for 10 6 Mq black hole 
radiating at 0.5Mem is ~ 60 eV. [The black hole mass was 
estimated by assuming a bolometric luminosity of ~ 10 x 
the 0.2-10 keV X-ray value and taking this to be the Ed- 



dington limit. This, together with the assumed accretion 
rate, is an acceptable value, since NLSls are thought to be 
low mass systems with high mass-accretion rates (Pounds, 
Done & Osborne 1995).] This is an upper limit to the ac- 
cretion disc temperature, since areas more distant from the 
black hole will be radiating at lower values. It can be seen 
that the BB of ~ 80 eV is broadly consistent with this value, 
while the 200 eV component could not be produced through 
thermal radiation from the disc. 

A more realistic model involves Comptonisation: soft 
photons from the accretion disc are up-scattered by hot, 
thermal electrons, possibly located in a corona above the 
disc. A two-temperature distribution leads to the formation 
of both the soft excess and the harder power law slope. To 
determine whether Comptonisation could be used to explain 
the Mrk 896 spectrum, the Xspec model comptt was used. 
Initially, just the soft excess was modelled in this fashion; 
i.e., the power law component was used for the higher energy 
part of the spectrum, with the iron line fixed as in Fit 2. This 
led to a very good fit (\l = 704/715), for T = 2.06 ± 0.04, 
together with an input photon temperature of (67 ± 5) eV 
and a Comptonised component of kT = (0.47 ± 0.23) keV 
and r = 6.7 ± 2.3. 

When using a second comptt component to replace 
the power law, it was found that the temperature for the 
hotter distribution could not be very well constrained. This 
is due to the fact that the exponential cut-off produced in 
the Comptonisation fit lies outside the XMM energy band- 
pass and cannot, therefore, be easily determined. However, 
a statistically good fit is obtained when the temperature is 
fixed at 200 keV; the model parameters are given as Fit 8 
in Table [5] and, again, the iron line was kept as in Fit 2. 
It must be noted that, because of fixing the hotter kT, the 
error on the optical depth appears small. 

There are many theories for the production of the X- 
ray spectra through Comptonisation (and others which do 
not invoke Comptonising distributions). When modelling 
the spectrum of Mrk 896 using COMPTT, it was assumed 
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Figure 3. An unfolded plot of the best fit thCompfe model; 
again, only the PN data are shown. The model consists of disc 
photons at ~ 60 eV Comptonised by two electron distributions, 
at ~ 0.3 and 200 keV respectively. 



that cool disc photons were Comptonised by one of two tem- 
perature distributions. Although this is a possible model, it 
is more likely that there is a 'layered' structure, such that 
most of the thermal photons are first Comptonised by the 
cooler of the electron distributions, forming the photons we 
see as the broad soft excess; following this interaction, some 
of these photons will be further Comptonised by the hotter 
electrons (likely to be produced through magnetic reconnec- 
tion above the accretion disc; this electron distribution may 
be non-thermal) , producing the power law tail seen at higher 
energies. If this is indeed the case, then the input photon 
temperature for the hotter Comptonisation should be that 
produced by the first Comptonisation stage. This is the situ- 
ation assumed when fitting with thCompfe, an alternative 
Comptonisation model (Zycki, Done & Smith 1999) which 
takes into account the roll-off of the power law at lower en- 
ergies due to the input photon distribution. In this respect, 
it is a more self-consistent model, and tends to give better- 
constrained results. Using thCompfe produces an excellent 
fit, with xl of 704/715 (Fit 9 in Tabled . The temperature of 
the hotter distribution was fixed to 200 keV, but the cooler 
electron temperature was allowed to float. Figure [3] shows 
the spectrum. It should be noted that assuming the same 
input photon population for both Comptonising distribu- 
tions also produces as good a fit, although the temperature 
of the electrons producing the soft excess photons is not well 
constrained. Overall the Comptonisation models provide as 
good a fit to the data as the BB model since they repro- 
duce the spectral curvature better than multiple power law 
models. Using the F-test, a value of F > 18 is obtained from 
the Comptonisation model over the power law fits of the soft 
excess; this corresponds to an improvement of > 99 per cent. 

3.1.3 Ionised disc fit 

An ionised accretion disc reflection model, described by Bal- 
lantyne et al. 2001, was also used, initially to fit the 2-10 keV 
energy range. The best fit obtained was found to consist of 
a power law component, of V = 2.04 ± 0.04, with an ionisa- 
tion parameter of 1.37 ± 0.35 erg cm s _1 (i.e., only weakly 



ionised) and a reflection component, R = 0.52 i ig, where 
R = 0,/2-k. This gave a reduced \ 2 value of 250/266, which 
is not significantly different than the value of 246/263 for a 
Gaussian fit. 

Ionised reflection may contribute to the soft X-ray cur- 
vature observed in broad-band AGN X-ray spectra, often re- 
moving the need for multiple blackbodies (e.g., for Mrk 205, 
Reeves et al. 2001; Mrk 359, O'Brien et al. 2001; Mrk 509, 
Pounds et al. 2001). Even upon including the ionised reflec- 
tion model in the fits, and allowing the ionisation parameter 
to vary, two BB components (with temperatures very similar 
to the fit without the ionised reflector) were still required to 
fit the broad-band Mrk 896 spectrum (xl = 716/717); how- 
ever, the normalisations of the blackbody components were 
reduced by approximately 25 per cent in strength. It is clear, 
therefore, that not all of the soft excess in this object can be 
explained purely through the reprocessing of the primary, 
hard X-ray component by an ionised disc. 

3.2 RGS data 

The RGS instrument allows a more detailed investigation 
of the soft X-ray spectrum. In an earlier paper (O'Brien et 
al. 2001), the NLS1 Mrk 359 was analysed. The RGS data 
for this object were found to show an absorption trough, 
possibly corresponding to Fe M ions (e.g., Sako et al. 2001), 
together with emission lines from O viii Lya and the Ne IX 
and O vil triplets. The RGS spectrum of Mrk 896 shows a 
weak absorption trough around 16.9 A (rest frame), which 
could correspond to an Fe M feature, and a small peak 
around the energy expected for O vm Lya (rest frame 
wavelength of 18.9 A). However, fitting a similar model to 
that used for Mrk 359 provides only equivalent-width up- 
per limits for these spectral features of EW < 18.4 eV and 
EW < 12 eV respectively. Similarly, the respective com- 
bined upper limits for the Ne IX (A res t = 13.6 A) and O vil 
(Aresi = 22.0 A) triplets are EW < 11.6 eV and <10.3 eV. 
There are also no strong absorption edges within the data, 
with t(0 vii) < 0.35 and r(0 vm) < 0.16 (at the 90 per cent 
confidence level). 

The lack of strong spectral features in the soft energy 
region implies that the observed soft excess is not dom- 
inated by a blend of soft X-ray lines (e.g., Turner et al. 
1991). The shape of the spectrum is also very different from 
MCG— 6— 30— 15 and Mrk 766, where relativisitc lines were 
used to fit the soft excess by Branduardi-Raymont et al. 
(2001). 

3.3 OM data 

Mrk 896 was observed using the V-band filter of the Optical 
Monitor only. A magnitude of my = 15.28 ± 0.01 was ob- 
tained, corresponding to a flux of 2.735 mjy. The OM uses 
a small aperture size of six arcsec. Winkler (1997), however, 
used a twenty arcsec aperture, which includes more of the 
host galaxy, so obtaining a brighter magnitude of ~ 14.4. 
Figure^Jshows the position of the optical points (from Win- 
kler 1997 and the OM - circle and star markers respectively) 
in relation to the X-ray data, together with the radio mea- 
surement obtained from Bicay et al. (1995) and four IR data 
points, obtained from the IRAS Faint Source Catalogue, ver- 
sion 2.0 (Moshir et al. 1990). 
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Table 2. 



Comptonisation fits to the broad-band XMM data. 
t tied to cooler component * frozen 



COOLER COMPTONISING COMPONENT 



HOTTER COMPTONISING COMPONENT 



Fit Model 



input photons Compt. temp. optical depth input photons Compt. temp. optical depth x 2 /dof 
(eV) (keV) (eV) (keV) 



8 comptt 66 ± 5 

9 thCompfe 58 ± 5 



0.54 ± 0.34 
0.28 ± 0.18 



6.2 ± 3.1 
13.8 ± 6.3 



66* 
280' 



200 / 
200- f 



0.14 ± 0.01 
0.77 ± 0.03 



704/715 
704/715 




14 

log v (Hz) 

Figure 4. A Spectral Energy Distribution of Mrk 896, showing 
the radio (triangle), IR (squares), optical (circles and star) and 
X-ray (crosses) measurements. 



4 DISCUSSION 

A simple analysis of the broad-band XMM-Newton spec- 
trum of Mrk 896 reveals a soft excess lying above a power 
law, fitted over the 2-10 keV band, as is typical for obser- 
vations of AGN (e.g., Pounds & Reeves 2002). The lack of 
strong spectral features in the RGS spectrum appears to rule 
out the possibility that the soft excess could be due to the 
blend of soft X-ray lines. It is found that a two-temperature 
Comptonisation model, where photons from the accretion 
disc undergo inverse-Compton scattering with hot electrons 
to produce the soft excess and hard 'power law' respectively, 
fits the spectrum very well. The input photons are ~ 60 eV 
(appropriate for emission from the inner-most regions of 
the accretion disc) and the Comptonising distributions have 
temperatures of kT ~ 0.3 and 200 keV. It is possible that 
the input photons for the hotter electron distribution are 
those which have been previously Comptonised by the soft- 
excess-producing population; however, the geometry of the 
inner regions of active galactic nuclei is unknown. It could 
be that photons emitted from the accretion disc are directly 
Comptonised to form the harder power law slope. Alterna- 
tively, the hotter electrons may be a non-thermal distribu- 
tion (see Coppi 1999 and Vaughan et ai. 2002 for details of 
hybrid thcrmal/non-thermal plasmas). 

Considering the temperatures found for the blackbody 
fit, the cooler of the two components (~ 85 eV) is consis- 
tent with being the high-energy tail of the big blue bump, 
while the hotter, ~ 225 eV BB corresponds to the tem- 



perature of the soft-excess-producing Comptonising com- 
ponent. (The power law replaces the higher-temperature 
Comptonised term.) 

It should be noted that it is not possible statistically 
to differentiate between the blackbody and Comptonisation 
models fitted to these data. 

A neutral Fe Ka emission line is also found. Although 
an ionised accretion disc model can be invoked to explain the 
slight broadening of the iron emission line, this does not ac- 
count for the entire curvature at lower energies; that is, the 
soft excess cannot arise simply through the reprocessing of 
the hard X-ray continuum. The line is, however, consistent 
with being a narrow feature; it is becoming increasingly ap- 
parent that almost all broad-line AGN, below a luminosity 
of ~ 10 45 erg s _1 , show this feature, with the lower lumi- 
nosity objects having the stronger lines. It is thought that 
this narrow line is formed through fluorescence in distant, 
cool matter, possibly the BLR or putative molecular torus 
(e.g. Ghisellini, Haardt & Matt 1994), with higher luminos- 
ity objects having smaller covering factors (O'Brien et ai. 
2001; Pounds & Reeves 2002). 
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